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Abstract We present the results of the INTEGRAL monitoring campaign on the accreting low mass X-ray binary 
pulsar GX 1+4 performed during the Galactic plane scan of the INTEGRAL Core Programme. The source 
was observed in different luminosity states ranging from I/20-40kcV = 1-7 x 10 _10 erg cm -2 s _1 , to I/20-40keV = 
10.5 x 10~~ 10 ergcm -2 s" 1 for about 779 ks from March 2003 until October 2004. Our observations confirm the 
secular spin down of GX 1+4 with the spin period (P s ) varying from 139.63 s to 141.56 s. In the highest luminosity 
state, a spin-up phase is observed. The phase-averaged spectrum of the source was modelled either with an 
absorbed cut-off power law or with a Comptonization model with significantly different parameters in the two 
brightest luminosity states. No evidence of any absorption-like feature is observed in the phase averaged spectrum 
up to 110 keV. At highest luminosity, the source is found to pulsate up to 130 keV. Phase resolved spectroscopy 
reveals a phase-dependent continuum and marginal evidence for an absorption feature at 34 ± 2keV in the 
descending part of the pulse. If interpreted as due to electron resonant cyclotron scattering, the magnetic field 
in the emitting region would be (2.9 ± 0.2) x 10 12 (1 + z) G where z is the gravitational red shift of the emitting 
region. We also observed a very low luminosity state, typical of this source, which lasted for about two days during 
which the source spectrum was modelled by a simple power law, and a pulsed signal was still detectable in the 
15-100 keV energy range. 
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1. Introduction 

The X-ray binary pulsars (XRBPs) were discovered 
more than 30 years ago wi t h the pioneer observation of 
Cen X-3 by iGiacconi et al Even though the ba- 

sic mec hanisms of such a pulsed emission were understoo d 
quickly <)Pringle fe Reesll972tlDavidson fe Ostrikeill973|) . 
XRBPs still present many puzzling aspects. In particular, 
their wide-band spectral behaviour has not been explained 
yet on the basis of a unique coherent physical model. 
Radiation from XRBPs originates from the accretion of 
ionized gas, from a nearby companion, an O or B star 
(M > 5M ) in High Mass X-ray Binaries (HMXRBs), a 
later than type A star (M < M©) or a white dwarf in 
low mass X-ray binaries (LMXRBs), onto the magnetic 
pole regions of strongly magnetized (B ^10 12 G) rotat- 
ing neutron stars. At several hundred neutron star radii, 



the plasma is threaded and then funneled along the mag- 
netic fiel d lines onto the magnetic poles at the neutron sta r 
surface l|Pringle fc Reesll972l:lDavidson fc Ostrikerll973T) . 
Pulsation is generated because of the nonalignment be- 
tween the magnetic and rotational axes. Spin periods are 
not constant and spin-up and spin -down behaviours ar e 
observed on a variety of time scales ( Bildsten et al.lll99^l . 
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According to accret ion theory l|Ghosh et alJ Il977t 
lOhosh fc Lambl |l979alhl ). the accretion process onto a 
neutron star surrounded by a disk strongly depends on the 
structure of the transition zone where the Keplerian disk 
meets the pulsar magnetosphere. The occurrence of accre- 
tion is mainly governed by the ratio between the magneto- 
spheric radius r m , where the ram pressure of a spherically 
symmetric inflow equals the magnetic pressure, and the 
corotation radius r co , where the magnetic field lines move 
at the local Keplerian velocity. 
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When r m < r co , the plasma is forced to stream along 
the field lines and eventually falls onto the neutron star's 
magnetic poles. When r m > r co , the accretion stops due to 
the onset of the centrifugal barrier, wh i ch is k nown as the 
prope ller effect Ijlllarionov fc Sunvaevl Il975l IStella et all 
Il986l) . If r m ~ r co , the system is in a near equilibrium 
state. In this condition, spin-up and spin-down periods, 
due to the transfer of angular momentum from the ac- 
creted plasma to the NS, can occur when variations in the 
accretion rate (M) take place. When M is sufficiently high, 
the star experiences a spin-up torque; as M decreases, 
the spin-up torque falls and vanishes at a critical accre- 
tion rate. For a lower accretion rate, the star experiences 
a spin-down torque whose magnitude increases until M 
reaches th e minimum accretion ra te consistent with steady 
accretion llGhosh fc Lamblll979bh . 

Such torque reversal episodes characterise the peculiar 
spin history of the (P - 140 sec) LMXRB GX 1+4. The 
source, which appears to accrete from the wind of the 18 
Mag, M5III companion, V2116 Ophiuchi l|Davidsen et alJ 
Il976|) . exhibited in the 1970s the highest spin up-trend 
ever obse rved in any p ulsar with vjv: as high as ~2% 
per year l|Nagaselll989l) . The orbita l period of the sourc e 
was reported t o be ~ 304 d a ys by IPereira et alJ l|l999h . 
More recently, iHinkle et alJ l|2003|) reported a P or b ~ 
1042 days, based on infrared spectroscopic properties of 
V2116 Ophiuchi. 

No observations were made from 1980 to 1982; later 
EXOSAT observations in 1983 and 1984 failed to de- 
tect the source, indicating an X -ray flux decrease of a t 
least two orders of magnitude (jHall fe Davelaarl Il983j) . 
In 1986, GX 1+4 wa s observed during a balloon flight 



I Grccnhill et al 1 11 98^ and exhibited a spin period of 
111.8s, indicating that one or more spin-down episodes 
occurred between 1980 and 1986. However, this episode 
must have been localized in time, since a collection of 
different data since 1980 shows that GX 1+4 was again 
spinning up at an average ra te of 6.6 x 10~ 8 s/s sim ilar to 
the rate it had in the 1970s l|Greenhill et al.lll989() . 

The pulsar reappeared with a reversed period deriva- 
tive and low lum inosity in 1987 l|Makishima et alJll988l 
iMonv et al1ll99lj) : and during 1989-1991, it continued to 
spin down rapidly with a mean rate ab out half, in absolute 
value , of the previous spin-up rate I Chakrabartv et aD 
I1997|) . 

On an irregular basis, GX 1+4 presents spin-up 
episodes with the absolute intensity similar to the spin- 
down trend. To explain such episodes on the basis of some 
apparent anti-correlation between torque and luminosity 
in the 20-60 keV band seen by BATSE, seve ral authors 
( Chakrabartv eTaTI Il997t iNelson et al] Il997|) have sug- 
gested the formation of a retrograde disk implying a mag- 
netic field of B ~ 10 12 G on the surface of the neutron star. 
Instead, other authors i Beurle et al.lll984l iDotani et alJ 
Il989t IMonv et al.lll99lh inferred the highest dipole mag- 
netic field of any accreting pulsar from the torque reversal 
episodes: B = (1-3) x 10 13 G on the surface of the neutron 
star. 



Cyclotron resonant scattering features (CRSFs) pro- 
vide a tool for direct measurement of the magnetic- 
field strengths of accreting pulsars since they appear at 
E cyc = 11-6 B12 x (1 + z) _1 keV, where B 12 is the mag- 
netic field strength in units of 10 12 G and z is the gravita- 
tional red shift of th e emitting region i Coburn et al.l l2002: 
Idi Salvo et"aT1l2004h . No CRSFs have been detected until 
today in the spectrum of GX 1+4. This lack of detection 
is not surprising if the estimate based on the standard ac- 
cretion torque theory is correct. In fact, the inferred mag- 
netic field of - 10 13 G converts to E cyc « 100-300keV, 
which is beyond the detection limit of hard X-ray satellites 
like Beppo-SAX a nd RXTE for the flux level and sp ectral 
shape of GX 1+4 llNaik et al.ll2005t lGallowavll2000l) . 

Moreover, evidence has been found that the hard- 
ness of the XRBP spectrum is closely related to the field 
intensity. Ginga, Beppo-SAX, and RXTE observations 
of binary accreting pulsars have revealed a correlation 
between line energy and spe ctral hardness in the high- 
energy part of the spectrum llMa kishima fe MiharaH l992; 
Idal Fiume et al.ll2000l ICoburn et al.ll2002j) . Since the ob- 
served spectrum of GX 1+4 is one of the hardest found 
among all X-ray binary pulsars, this relation suggests a 
magnetic field in excess of 10 13 G. Similar re sults have 
been found for GS 1843-00 llPiraino et al.ll2000l) . 

In intermediate and high luminosity states 
(^2-iokcV = 10 -10 -10~ 9 ergcm -2 s" 1 ), its emission 
has been described either with a cut-off pow er law 
llKotani et al.l FflflS EM ICini fe Smithl 

iRea^^Ll 120051) or with a Comptonization mode l 
llGallowayll27)0(riGallowav et alJl200l( iNaik et allEoOfih . 
In both cases the source shows high intrinsic absorption 
and exhibits a strong Iron line at 6.5-7 keV with an 
equivalent width of 0.2-0.5 keV. 

On an irregular basis, GX 1+4 presents low luminos- 
ity states in which the flux decreases to a few L2-iokcV = 
10 -11 erg cm -2 s _1 . In these states, even though the other 
spectral characteristics are highly variable, the source al- 
ways shows an enhanced Fe emission at 6.5-7 keV (few 
keV equivalent width) and a harder spectrum in which 
the cut-off often cannot be found or is at higher en- 
ergy with respect to the b r ight state faea et alJ l2005t 
INaik et all Eoosl ICul Il997t ICui fc Smithl l2004h. These 
episodes can last from about one da y l|Gallowav et alJ 
Eoorl lRea et al]l2005UNaik et a,lJl2005h to several months 
l|Cuilll997l: ICui fe Smithll2004tl. The phenomenol ogical in- 
terpretation is ambiguous: ICui fe Smithl l)2004() explain 
the low luminosity episodes as due to the propeller effect, 
while the residu al power law emission i s produced by accel- 
erated particles : lGallowav et alJ 1 2000i) p resent a model in- 
volving the accretion column disruption; iRea et all l|2005h 
suggest that the pulsar is hidden by a thicker accretion 
disk. In the last case, the residual radiation is produced 
by the reflection of the disk edge and by partial trans- 
parency of the medium at high energies (> lOkeV). 

The main goal of this paper is to study the spec- 
tral and timing characteristics of GX 1+4 observed 
over a wide range of luminosity, within the INTEGRAL 
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Galactic Plane Scan (GPS) part of the Core Programme 
^Winkler et all l2003|) from March 2003 until April 2004. 
The search for CRSF is also one of the key goals of our 
study. 

In Sect. |2 we describe the observations and the data 
analysis method. In Sect.|3|we report on the results ob- 
tained from analysing light curves, spin period, pulse pro- 
files and phase-resolved and phase-averaged spectra. In 
Sect. Q]we discuss our results in relation to previous ob- 
servations. Finally, in Sect. [5] we draw our conclusions. 

2. Observations and data analysis. 

The European Space Agency's International Gamma- 
Ray Astronomy Laboratory (INTEGRAL), launched in 
October 2002, carries 3 co-aligned coded mask telescopes. 

- The imag er IBIS (Imager on B oard the INTEGRAL 
Satellite; lUbertini et all 12003ft . which allows for 
12arcmin angular-resolution imaging in the energy 
range from 15keV to 600 keV with energy resolution 
of 6-7%; IBIS is made of a low-energy (15 -600 keV) 
CdTe detector ISG RI llLebrun et al.l l2003h and of a 
Csl layer, PICsIT l|Labanti et al.H2003ft . designed for 
optimal performance at 511 keV, with an overall nom- 
inal energy range between 175 keV and lOMeV. 

- The spectrometer SP I (Spectrometer on INTEGRAL; 
IVedrenne et~al"ll2003h . sensitive from 20keV to 8MeV 
with an angular resolution of 2.5 deg and an energy 
resolution of a few times 10 -3 . 

- The X-r ay monitor JEM- X (Joint European X-ray 
Monitor; iLund et al.l l2003) made of two independent 
units JEM-X1 and JEM-X2, sensitive from 3keV to 
34 keV, with an angular resolution of 1 arcmin and an 
energy resolution of 10-15%. 

As part of its core programme, INTEGRAL performs 
regular scans of the Galactic plane (GPS), with the aim 
of monitoring the timing and spectral properties of the 
known X-ray sources, discovering new transient sources, 
and mapping the di ffuse emission of the Galactic plane 
llWinkler et al.ll2003l) . 

From March 21, 2003 to September 17, 2004 the region 
containing GX 1+4 was observed several times. In order 
to have the shadowgram of GX 1+4 fully coded on the 
ISGRI detector plane, only pointings where GX 1+4 is 
observed at less than 6 deg off-axis were used. This gives 
the best S/N performance for ISGRI. 

The data analysed in this paper consist of 435 sci- 
ence windows for a total observing time of 779 ks and 
are divided into four intervals of semi-contiguous science 
windows as shown in the light curve of Fig. We de- 
fine time intervals as follows: MJD 52 720-52 750 observa- 
tional period 1 (OPl), MJD 52915-52925 (OP2), MJD 
53050-53090 (OP3), and MJD 53 250-53290 (OP4). We 
also distinguish one sub-period corresponding to a low- 
flux episode in the middle of an intermediate flux state 
(OP3L), which started on MJD 53 067.25 and ended on 
MJD 53 068.75. In OPl, OP2, and OP3L, the source flux 
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Figure 1. Light curve of GX 1+4. We plot the 
INTEGRAL ISGRI count rates in the 22-110keV band 
during the galactic plane scan pointings with an off-set an- 
gle less than 6 deg. The binning is one science window cor- 
responding to 1.6-3.6ks. The time scale is modified Julian 
days (MJD=JD-2 400 000.5); MJD 52700 is 1 March 2003. 



was too low for a detection with JEM-X or SPI. In OP4 
and in the last part of OP3 (after MJD 53 072), it was pos- 
sible to use only JEM-X1, while in the first part of OP3 
only JEM-X2 was available. In OP3 we did not analyse 
the 14 ks of JEM-X1 data due to the low statistics; the 
reported exposure refers to JEM-X2. The lower exposure 
time of JEM-X compared to ISGRI is due to the smaller 
field of view of this instrument. The complete list of the 
observations is reported in Tabled 

We used version 5.1 of the INTEGRAL off-line stan- 
dard analysis software (OSA 5.1), officially distributed by 
the INTEGRAL Science Data Centre (ISDC), to obtain 
the phase-averaged spectra of SPI and JEM-X. For the 
latter instrument we limited the analysis to the energy 
range 4-21 keV, for which the JEM-X team ensures a re- 
liable calibration as stated in the OSA handbook. All the 
information regarding the official software and the stan- 
dard calibration can be found on the ISDC web-site at 
http : //isdc . unige . ch/. 

To analyse the ISGRI data we had to face two kinds 
of problems that are not optimally solved by the OSA 
pipeline: the lack of tools for performing straightfor- 
ward timing and phase-resolved analysis and a not-yet- 
optimized energy reconstruction of the events. This re- 
quired the development of specific tools, which are de- 
scribed in the appendix. 

3. Results 

3.1. Light curves 

The ISGRI light curve of GX 1+4 in the 22- 
110 keV energy range, averaged over the science win- 
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Table 1. Summary of the INTEGRAL observations of GX 1+4. In the table we show the name for the observational 
period, its start and end time, the exposures of the three INTEGRAL instruments used here and the corresponding 
total number of source counts. 



name 


start [MJD] 


end [MJD] 


exposure [ks] 






Source counts 










ISGRI 


JEM-X 


SPI 


ISGRI 


JEM-X 


SPI 


OP1 


52 720.8076 


52 748.9027 


72 






2.9 x 10 b 






OP2 


52 917.6381 


52 921.7215 


97 






5.2 x 10 5 






OP3 


53 052.0923 


53 085.5729 


155 


47 


137 


1.9 x 10 6 


2.0 x 10 5 


4.2 x 10 3 


OP4 


53 252.4430 


53 287.3416 


423 


154 


351 


2.4 x 10 6 


3.8 x 10 5 


3.9 x 10 3 


OP3L 


53 067.2531 


53 068.7772 


32 






9 x 10 4 
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Figure 2. Light curve of GX 1+4 observed by ISGRI dur- 
ing the GPS for OP3 in the 22-110 keV energy band. 
The region with the dashed border highlights the low- 
luminosity episode, which is reproduced in the inset in 
the upper right corner. The binning is one science win- 
dow corresponding to about 1.7-3.6 ks depending on the 
operational mode. Errors are at la level. 



dow exposure (typically 1.6-3.6ks), is shown in Fig. ^ 
During the INTEGRAL AOl Core Program Monitoring, 
GX 1+4 evolved from a weak state (£20-40 kcV — 2 x 



! s 1 ) to a much brighter one (L 



20-40 kcV 



10 10 ergcm 

10 x lO-^ergcm^s- 1 ). 

During OP3, from February 17, 2004 to March 23, 
2004, the source exhibited erratic variations of about one 
order of magnitude in flux on a few kilo-seconds time 
scale (see FigEJ. In particular, from MJD 53 067.25 until 
MJD 53 068.75, GX 1+4 underwent a very low luminosity 
state (OP3L) observed by INTEGRAL for -32 ks. 

Since these low-lumin osity episodes are peculiar in the 
behaviour of the source llRea et a,lJl20fii lGa,llowa,vl l2000: 
ICuifc Smith1l2004h . OP3L data are analysed separately 
and excluded from the average analysis of OP3. 



3.2. Spin period 

Spin frequencies in the four observational periods were 
determined separately. First, we computed a preliminary 
power spectrum from the light curve binned with a 10 s 
time resolution. The ~ 0.007 Hz peak position gave us a 
tentative period from which we produced a pulse profile 
Pt{4>) for each science window. Then, we computed the 
phase of the maximum of the pulse first Fourier compo- 
nent in each science window: 



0o (t) =tzn- 1 (I c {t)/I s {t)) , 



where 



I c (t)= / P t ((A)cos(0)d0 and 



/,(*) = / P t {4>) sin($(ty 



(1) 
(2) 
(3) 



and t is the science window central time. 1 If the law 
we used for folding is not optimal, <f>o(t) is not con- 
stant throughout the observational period; however, we 
can derive the optimal folding law from its variation. Let 
Ac/)f(t) = 4>o(t) — <fio(to) be the phase shift measured after 
folding the data with the temptative period; we recon- 
struct the actual phase shift A</>(i) by fitting A</>(i) with 
a second-order polynomial 



A0(t) = Q + a(t-t o ) + b(t-t Q ) 2 



(4) 



where 4>q is the starting phase shift determined by the fit 
and to the reference folding time chosen close to the begin- 
ning of the observational period. The coefficient a is used 
to correct the tentative choice of the folding frequency v t 
as 

u c = u t - a, (5) 

where v c is the corrected pulse frequency. The coefficient 
b is linked to the pulse frequency derivative during the 
observational period by the relation 



= -26. 



(6) 



1 For the pulse shape of GX 1+4, most of the power is con- 
tained in the first harmonic; for sources with a double peak 
this method should be used with caution. 
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Figure 5. The background subtracted pulse profiles for 
OP2, OP3, and OP4 in the energy ranges indicated in 
the plot. The reference time, spin period, and spin period 
derivative are listed in Table [21 



Finally, we verified that A0(t) obtained with the corrected 
folding law is constant within the errors and therefore the 
derived parameters are correct. A visual inspection of the 
residuals in Fig. [31 confirms that it is not necessary to 
use a higher-order polynomial and that the connection of 
phase shifts separated by wide gaps within the OPs is 
self consistent, even though in principle it could be non- 
unique. The different OPs cannot be phase-connected due 
to their time separation and to the complex spin evolution 
of the source, with e.g. a torque reversal between OP3 and 
OP4. 

Due to the low number of source counts for OP1, spin 
periods and their derivatives were only measured for OP2, 
OP3, and OP4. In Table[21we summarise the values for the 
spin period and its first derivative, and the reference time 
at which they are measured. The relative time intervals are 
specified in Tabled In some science windows in OP2 and 
OP3, the source was too weak for a good determination 
of the phase of the maximum, we so excluded them from 
the analysis, which produced an exposure time reduction 
of less than 5%. The source shows the longest spin period 
of its observational hist ory P s —141 s, which con firms the 
overall spin down trend IjOhakrabartv et alll 997). In OP4 
we find a sign inversion of the spin period derivative, which 
is the signature of a torque reversal episode. 

The secular spin-down trend we measured from a 
linear fit to the three measured pulse period values is 
P = (6.6 ± 0.8) x 10~ 8 s/s, but from the analysis of phase 
shifts we find local spin period derivatives that are larger 
in absolute value. This is consis tent with the source be- 
haviour in previous observations ( Chakrabartv et alJl997t 
IPereira et alJl999^ : local spin-period variations are super- 
imposed on a global spin down trend. 



3.3. Pulse profiles 

Using the values of P and P in Table [21 we extracted the 
background subtracted pulse profiles of ISGRI for differ- 
ent energy intervals using the tools described in the ap- 
pendix B. The pulse profiles for OP4 are shown in Fig. 0] 
the ones for OP2, OP3, and OP4 in Fig. Phase zero 
was arbitrarily chosen at the bottom of the dip in the 
15-25 keV profile for OP4 and at the pulse minimum for 
the other OPs. At the highest luminosity, the source is 
observed to pulsate up to 130 keV. 

As a first approximation all pulse profiles show a sin- 
gle broad peak - as stated repeatedly in the literature. 
There are, however, clear substructures. A dip like struc- 
ture within the broad minimum (around pulse phase zero) 
is evident at energies up to 60keV (Fig.|3}. This is a well- 
known featur e already found in RXTE and Beppo -SAX 
data (see e.g. lGallowav et allEoOll iNaik et alJl2005h . The 
pulse appears asymmetric and looks like a superposition 
of two (or more) peaks. This is also strongly suggested 
by various pulse profiles found in th e literature (see e.g. 
iGiles et allbood ICui fc SmithlbOO^ . Following this sug- 
gestion, we modelled our pulse profiles of OP4 by the sum 
of two Gaussian components (leaving out four data points 
that define the dip around phase 0). The model fits are 
shown as solid lines for all profiles below 100 keV in Fig.QJ 
Each Gaussian is characterised by three parameters: cen- 
troid, width, and amplitude (a constant flux was set to 
zero since it was found to be unconstrained when it was 
left as a free parameter) . Most fits are quite good (reduced 
chi-squared — 1), and others are acceptable (reduced chi- 
squared up to 1.4 for 40 d.o.f.). Our findings from these 
fits for OP4 (Fig. 0} are the following: 

1. The centroids of both peaks are consistent with a con- 
stant in phase throughout the entire energy range: 0.39 
and 0.71, for the first and second peaks, respectively. 

2. The widths of both peaks are slightly variable: while 
for energies above 40 keV both peaks have a width con- 
sistent with 0.18 phase units, at lower energies the first 
peak becomes narrower while the second peak becomes 
wider (reaching 0.16 and 0.26 phase units, respectively, 
for 15-25 keV). 

3. As evident from the visual inspection and from the 
ratio of the amplitudes of both peaks, the second peak 
is the dominating one. The ratio of amplitudes (first 
peak to second peak) increases from 0.5 at 15-20 keV 
to 0.9 at 40-50 keV with a subsequent small decrease 
to -0.7 at 80-100 keV. 

4. There is only a weak indication of a third small peak 
inside the pulse minimum around phase 0.1 (see, for 
instance, the 30-35 keV profile), reminisc ent of the 
"bump" detected bv lGreenhill et alJ l)1989f) and there- 
after confirmed in othe r observations (e.g. iNaik et ail 
l2005tlGiles et aLlEoOoh . 

Pulse profiles of OP2 and OP3 (Fig. 0) confirm the 
general findings from OP4, although in OP2 and OP3 the 
two subpeaks contribute almost equally to the main pulse 
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Figure 3. The phase shifts are obtained using a constant folding period with the best-fit second-order polynomial 
overplotted (upper panels) and the residual phase shifts obtained with the folding parameters of Tabled (lower panels). 
The phase shifts can be converted into time units by multiplying by the folding period of 140 s. Left is for OP3 and 
right for OP4. The time scale is in seconds after the reference times of Table [21 



Table 2. Values for the spin period (P) and spin period derivative (P), in the three observational intervals where 
the source was bright enough for this determination, are reported. The errors are at la level computed from the fit to 
the phase shifts. 



conv. name 


reference time [MJD] 


P [s] 


P[s/s] 


OP2 


52917.638038 


139.630 ± 0.006 


(1.5 ±0.1) x 10"' 


OP3 


53052.091969 


140.6132 ± 0.0002 


(1.034 ±0.002) x 10 -7 


OP4 


53252.442934 


141.56488 ± 0.00014 


(-0.5331 ± 0.0004) x 10~ 7 



at all energies. An impressive demonstration of the exis- 
tence of different emission components in G X 1+4 and 
their variability is given in Figs. 4 and 5 of Idles et alJ 
( 2000) where three peaks (and the dip structure) are very 
apparent. 

To obtain the pulsed fraction we fitted the pulse pro- 
files with a sine function, plus a constant and computed, 
from the sine amplitude A and the constant C, the pulsed 
fraction as 2A/ (A+C). Pulsed fractions in different energy 
intervals (where the S/N permitted such an analysis) are 
reported in Fig.^l Interestingly, the intermediate luminos- 
ity observation (OP3) shows an average higher pulsed frac- 



tion in comparison with the brightest state (OP4). Pulsed 
fractions observed in OP2 are consistent with the others 
due to the large errors in the determination. 



Using the spin period found for OP3, we extracted 
a pulse profile for the low-luminosity episode (OP3L) in 
the 15-100 keV energy band. From the pulse profile of the 
right panel of Fig.[7]we can infer a pulsed fraction of 0.59± 
0.14 using the method outlined above. 
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Figure 4. The background subtracted pulse profiles of OP4 with energy ranges indicated in the figure. The vertical 
scale is in counts per second, the reference time, spin period, and spin-period derivative are listed in Table[5] The solid 
lines are the best-fit pulse profiles, the dotted and the dashed lines are the two Gaussian subcomponents used in the 
fit. 



3.4. Spectral analysis Broad-band spectral analysis, combining data from 

JEM-X, ISGRI, and SPI, was feasible only for OP3 and 
Spectral analysis was performed separately in the four OP4, since the source was not detected by SPI and JEM- 
OPs. Results for phase-averaged spectra are reported here x in the lower luminosity states. Thus, we firstly analysed 
while in Sect.ESwe describe our results on phase-resolved the ISGRI data alone for all luminosity states, and then 
spectroscopy. 
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Figure 7. The ISGRI spectrum (left) and the pulse profile (right) in the low-luminosity episode OP3L. The spectrum 
extends from 17 to 150 keV with no break and is fitted by a power law with a slope of 2.5 ± 0.2. The pulse profile is 
obtained in the 15-100 keV energy band. 
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Figure 6. The pulsed fraction in the three time intervals 
considered (see Table EJ. The crosses refer to OP2 (weak), 
the squares to OP3 (intermediate) , and the circles to OP4 
(high luminosity); the errors are la. 



we extended our analysis to the other instruments for the 
higher luminosity states. 

The ISGRI spectra were fitted in the 17-110keV 
energy range using a cut-off power law (N(E) oc 
E~ T cxp[— E/E c \) with the best-fit parameters sum- 
marised in Table |3J In OP1 and OP2 the parameter de- 
termination is poor due to the weakness of the source, so 
the fit results can be considered as generally consistent 
with the others. Instead, we note a significant difference 
between the best-fit parameters of OP3 and OP4, which 
can be studied better by adding the data from JEM-X 
and SPI. 



Table 3. ISGRI spectral-fit parameters for the model 
E~ r exp(—E/E c ) in the considered observing periods; 
OP3L is the low-luminosity episode. The flux is computed 
from the ISGRI data in the 20-40 keV energy range and 
expressed in units of 10 _10 ergcm~ 2 s . The uncertainties 
are at 90% c.l. 





Flux 


r 


E c [keV] 


X 7d.o.f. 




(20-40 keV) 








OP1 


1.7 ±0.5 


0.8 ± 1.0 


22 + f 


66/66 


OP2 


2.2 ±0.4 


0.9 ±0.6 


29+1 6 


55/66 


OP3 


4.9 ±0.3 


0.97 ±0.17 


27±* 


104/66 


OP3L 


1.8 ±0.5 


2.5 ±0.2 


> 150 


90/74 


OP4 


10.5 ± 0.2 


0.34 ±0.05 


20.1 ± 0.6 


97/66 



The broad-band JEM-X, ISGRI, and SPI (4-110 keV) 
spectra for OP3 and OP4 were fitted first using a cut-off 
power law with low-energy absorption ( WABS model in 
XSPEC) plus a Gaussian emission line at ~ 6.5 keV to 
model the blended neutral and ionized iron K-emission. 
The Gaussian emission line has parameters given by: 

N g / J2nag exp(—(E — £ , g ) 2 /4cTg). Best-fit parameters are 



shown in Table 0] (left panel) . The spectra were also suc- 
cessfully fitted with a Comptonization model (compTT 
in XSPEC) based on a sp herical geometry using the an- 
alytical approximation of iTitarchukl (^994) plus an ad- 
ditive Gaussian emission-line for the ionized K-emission 
from Iron and low-energy absorption. Best-fit parameters 
for this model are shown in Tabled (right panel). The 
compTT model was used with approximation 2 and for 
red-shift 0. kT is the temperature of the seed photons for 
the Comptonization, kT e the temperature of the scatter- 
ing electron cloud, while r is its optical depth in spherical 
geometry. The cut-off model is slightly better for OP3, the 
compTT model may be preferred for OP4. 
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Figure 8. Count-rate spectra (upper panels) and energy-flux unfolded spectra (lower panels) for OP3 (left) and OP4 
(right) of GX 1+4. The model for the continuum is compTT plus a low-energy absorption and Gaussian emission line 
at ~ 7keV. The data between 4 and 21keV are from JEM-X, and the data extending from 19 to 100 keV are from 
ISGRI. 



Significant differences between the spectra of OP3 and 
OP4 (at more than 90% c.l.) are found in the power law 
index and cut-off energy (cut-off model) and in the elec- 
tron temperature and optical depth (compTT model). In 
Fig.|Hlwe plot the broad-band count-rate and energy-flux 
unfolded spectra for OP3 and OP4 using the compTT 
model. 

The visual inspection of the fit residuals (see Fig.[S]up- 
per right panel) shows low points around 34 keV that are 
reminiscent of a similar structure at a 2.5a le vel seen at 
the s ame energies in Beppo-SAX observations l|Rea et alJ 
2005). However, the introduction of a multiplicative 
Gaussian absorption line -/V a b s exp(— (E— E^) 2 /2a 2 ) does 
not improve the fit, since its significance is much less than 
3a. 

The ISGRI spectrum of the low-luminosity episode 
(OP3L in Table was also extracted, but JEM-X and 
SPI did not detect the source in this state. The spec- 
trum (left panel of Fig. CJ was modelled by a power law 
with a slope of 2 5 ± 0.2 giving y 2 /d.o.f.= 9 0/74. The 
PEXRAV model ijMagdziarz fc Zdziarskilll995h proposed 



by iRea et all l|2005|) , which includes a reflection compo- 
nent besides the power law, is not constrained by our data. 

3.5. Phase resolved spectroscopy 

The OSA software (up to the current 5.1 version) does not 
provide tools for performing phase resolved spectroscopy 
for ISGRI and SPI. In the case of JEM-X the task is 
very cumbersome. Thus we developed dedicated tools (de- 
scribed in the Appendix lB| to analyse the ISGRI data. 

Phase-resolved spectra of OP3 and OP4 were extracted 
according to the phase interval selection shown in the up- 
per panel of Fig. 1101 and modelled with a cut-off power law. 
In OP1 and OP2 the source was too weak for this kind of 
analysis. As far as OP3 is concerned, best-fit parameters 
do not change with the pulse phase within the statistical 
uncertainties. On the contrary, the best-fit parameters of 
phase-resolved spectra of OP4 (Table change signifi- 
cantly with pulse phase: the power-law slope is not con- 
stant at a 99.9% confidence level, while the cut-off energy 
is not constant at a 90% confidence level. 
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Table 4. Results for the simultaneous fit of JEM-X, SPI and ISGRI spectra. The exposures are reported in Tabled 
and the uncertainties of the best-fit parameters are at 90% confidence level. The physical units of the parameters are 
reported between squared brackets except for the unabsorbed flux, which is computed in the 2-60 keV energy range 
and expressed in units of 10 -9 ergcm -2 s" 1 . 



Cut-off p 


awer law model 




compTT model 




parameter 


OP3 


OP4 


parameter 


OP3 


OP4 


N H [10 22 cm' 2 ] 


10 ±3 


5.0 ± 1.5 


N H [10 22 cm" 2 ] 


< 10 


25 ±2 


E g [keV] 


6.8± ;I 


6.5 ±0.1 


E g [keV] 


6.5 ±0.3 


6.6 ±0.3 


a g [keV] 


i±L 


< 0.4 


o g [keV] 


0.5 ±0.5 


0.2 fixed 


N s [10~ 3 cts cm" 2 s" 1 




2.7 ±0.5 


N s [10 -3 cts cm -2 s" 1 ] 


4±2 


1.1 ±0.6 


E c [keV] 


32 ±4 


21.8 ±0.7 


feT [keV] 


1.7 ±0.3 


1.6 ±0.2 


r 


1.2 ±0.2 


0.51 ±0.05 


kT e [keV] 

T 


15.4 ± 1.4 
5.1 ±0.5 


13.1 ± 0.2 
6.80 ±0.15 


Flux 2-60 keV 


1.77 ±0.07 


2.32 ±0.03 


Flux 2-60 keV 


1 7n+ - 04 

' u -0.16 


2.16 ±0.04 


X 2 /d.o.f. 


212/212 


311/199 


xVd-o.f. 


224/211 


303/201 




20 50 100 

Energy [keV] 

Figure 9. In the upper panel we show the ISGRI energy- 
flux unfolded spectrum extracted in the phase interval 0.6- 
0.85 The best- fit model has a multiplicative Gaussian ab- 
sorption line at (34 ± 2) keV. The other panels show, from 
top to bottom, the residuals without the introduction of 
the line and the residuals with its introduction. 



Although not formally requested by the \ 2 01 the fit, 
we also searched for absorption-like features in the phase- 
resolved spectra. This was, again, motivated by the resid- 
uals observed in the spectrum of phase 0.6-0.85 (Fig. EJ, 
the descending part of the pulse, and by the wea k ev- 
idence of such a line reported by iRea et all l|2005[) and 
iNaik et all l)2005|) . Using the Gaussian absorption model 
(see Sect. 13. 41) . the line centroid is found at 34±2 keV. The 
line width is not constrained by the fit and is set to 4 keV, 
compatible with the width of several absor ption lines ob- 
serve d at similar energy in other sources l)Coburn et alJ 
120021) . The line significance using an F-test is 2.4er, but, 
due to the questionable reliab ility of this statistical test 
in astronomical data analysis ( Protassov et alJl2002jk we 



conclude that the evidence for such a line, if any, is weak. 
However, systematic effects can be excluded because such 
a feature is not present in other pulse phases. 
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Figure 10. The pulse profile in the 25-50 keV energy 
band with the phase selection (upper panel), the slope 
(middle panel), and the cut-off energy (bottom panel) of 
the phase-resolved spectra of ISGRI in OP4. The best-fit 
parameters are obtained for the cut-off power-law model 
E~ r exp(E / E c ) without the introduction of the Gaussian 
absorption line in the energy range 17-110keV. The uncer- 
tainties represent the joint confidence limit at 90% prob- 
ability (Ax 2 = 2.7). 



4. Discussion 

The main results of our analysis of the INTEGRAL ob- 
servations of GX 1+4 are: 
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Table 5. Fit results for phase-resolved spectra, A^ 2 is referred to the passage from 72 to 70 d.o.f with the introduction 
of a Gaussian absorption line. The a of the Gaussian absorption line is fixed at a value of 4keV. The line significance 
is obtained using the F-test probability value. The uncertainties are the joint confidence limit at 90% probability. 



phase 0.10-0.35 0.35-0.6 0.6-0.85 0.85-1.10 

r 0.34 + 0.14 0.10 + 0.08 0.4 + 0.1 0.6 + 0.2 

E c [keV] 19 + 1 18.1 + 0.7 22.5 + 1.5 20 + 3 

£ abs [keV] 34 + 2 

iVabs 1.11+0.05 

Ax 2 69(72) -> 56(70) 

line significance 2.4<t 



1. The measurement of the longest spin period in the ob- 
servational history of GX 1+4: P s ~141s; 

2. The measurement of the secular spin-down trend 
P = (6.6+0.8) x 10 -8 s/s, in agreement with previous 
findings; 

3. The detection of a variation from spin-down to spin- up 
associated with the high luminosity state; 

4. Weak evidence in the phase-dependent spectra of 
an absorption feature with fundamental energy at 
33 — 35 keV that could be due to cyclotron resonant 
scattering; 

5. The observation of a low- luminosity episode where the 
source is found pulsating and the spectrum does not 
show evidence of any reflection component. 

In our Integral analysis we have also found, in line with 
previous observations from other satellites, evidence of en- 
ergy dependent pulse profiles, spectral variations between 
luminosity states, and phase dependent spectra. 

We confirm that the source is steadily spinning-down 
and maintaining a much lower luminosity than in the 
1970s. The values of the measured spin-period deriva- 
tive fall in the same range obtained on the basis of 
BATSE data. However, we observed a spin-up of the 
source in the high luminosity state: we exclude the pos- 
sibility that the spin derivative inversion is due to a 
binary-sys tem-induced Dopple r shift, as already demon- 
strated bv IPereira et al.l (U999J); it is more likely a signa- 
ture of a tor que-reversal episode, whi ch is rather common 
in GX 1+4 ijChakrabartv et aTlll997|) . In our observation 
the torque-reversal episode happened at a luminosity l evel 
compatible to the ones that iGhakrabartv et afl l(l997h re- 
ported as occurring at MJD 49 300, 49 700, and 50 300. 
Unfortunately, on the basis of our INTEGRAL observa- 
tion, it is not possible to narrow down the period at which 
the torque reversal occurred. It also needs to be stressed 
that a torque reversal does not imply such a luminos- 
ity value and vice - versa as can be seen from Fig. 2 of 
IChakrabartv et al I lll997h . where one sees spin-up episodes 
at lower luminosity and strong luminosity increases with- 
out a corresponding spin-up. 

Considerations of the nature of the emission mecha- 
nism can be made by discussing the pulse-profile evolu- 
tion with energy, together with the phase resolved spec- 
troscopy. When the source is brighter (L 2 o-40kcV = 10.5 x 



10 _10 ergcm~ 2 s _1 , OP4), the pulse profile can be heuris- 
tically described by two Gaussian sub-peaks centred at 
phases ~ 0.4, ~ 0.7: the relative intensity between the 
first and the second peaks increases with energy (Fig. . 
At ~ 20keV the pulse appears more peaked at an early 
phase while at energy just below 80keV the comparison 
with Fig. El shows that the two sub-peaks have equal in- 
tensity and the pulse shape becomes similar to the lower 
luminosity states observed in OP2 and OP3. 

As the pulse shape changes with energy, the spectral 
properties evolve with phase: in the ascending and top 
part of the pulse the spectrum is well-fitted with a power 
law of slope ~ 0.2 and cut-off energy ~ 19keV, while in 
the descending part the spectrum softens. This suggests 
a larger optical depth for scattering in the first sub-peak 
with respect to the second sub-peak, which turns out to be 
more similar to the phase-averaged spectrum of OP3. Such 
a phase-dependent behaviour is different fro m the one 
observed in previous phase-resolved analyse s I Gallowavl 
1200(1 iGallowavet "all l200lt iNaik et aljEool . which are 
focused on an absorption enhancement in correspondence 
with the pulse dip (at phase zero in our pulses) and do 
not find significant changes in the continuum within the 
remaining part of the pulse. 

The broad-band JEM-X/ISGRI/SPI spectra arc mod- 
elled either with an absorbed cut-off power law or with a 
Comptonization model. In OP1 and OP2, the fit param- 
eters are not well-constrained. In OP4 (when the source 
is brighter), we find a harder emission up to ~ 40keV, 
followed by a faster decay, in comparison to OP3. In con- 
trast to our results, ICui fc SmithI l(2004l) find that when 
the source is bright (~ 10~ 9 erg cm -2 s _1 in the 2-60 keV 
energy band), the spectral index is very stable between 1 
and 1.3 with e-folding energy at 28-40 keV. They also find 
that the fit parameters vary erratically with no apparent 
correlation to the flux on a time scale of a few days when 
the source is dimmer (spectral index between 0.2 and 2 
and e-folding energy between 11 and 45keV). These con- 
siderations imply that the spectral properties of GX 1+4 
are not determined uniquely by its luminosity state, but 
also depend on some still unknown parameter. 

The picture is slightly clearer when we use a ther- 
mal Comptonization model: the value of the seed com- 
ponent temperature (A: To ~ 1.6 keV) is always similar 
to the one found with Beppo-SAX and RXTE, while 
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the optical depth and plasma temperature are about the 
same but cannot be straightforwardly compared, confirm- 
ing the variability of the accretion mechanism within a 
stable system. Limiting ourselves to our data, where the 
source is brighter, the optical depth is higher with re- 
spect to the lower luminosity state, while the tempera- 
ture of the Comptonization plasma is lower. This spectral 
evolution with luminosity is consistent with the thermal 
Comptonization being t he dominant emission process (see 
e.g. the conclusions in iBecker fc W olff 2005): using this 
model the larger accretion rate can explain the enhance- 
ments of both luminosity and optical thickness in the ac- 
creting plasma. 

We observe weak evidence (at 2.4<r level) of an absorp- 
tion feature at 34 ± 2 keV in the descending part of the 
pulse profile, at phase [0.6-0.85] (see Fig.^J. Interpreting 
the absorption as due to the cyclotron resonant scatter- 
ing of electrons, the magnetic field in the emitting region 
would be (2.9 ±0.2) x 10 12 (l + .z) G where z is the gravita- 
tional red shift of the emitting region. Its presence or en- 
hancement just in the descending part o f the pulse is typi- 
cal of other systems such a s 4U 0115+63 
19991 iHeindl et a.lJ Il999h or Cen X-3 



Santangelo et al 



Santangelo et al 



1998). The inferred magnetic field would be one order 
of magnitude less than expected from the standard ac- 
cretion torque theory, but consistent with the one in- 
ferred ft^m_tIie_jnodelof_emission due to a retrograde 
disk ijChakrabartv et al"1 IT997). Although the calibration 
uncertainties of ISGRI suggest that this weak evidence is 
an indication rather than a measurement, we note that 
a similar indication comes from PDS on board of Beppo- 
SAX. On the other hand, it is not surprising that RXTE 
observations did not show spectral structures, although 
observati ons were carried out many times and i n a variety 



of states 



The sin- 



gle observations were shorter than those by INTEGRAL, 
and it is not possible to add many observations together 
because of the spectral variability of the source. 

The low-luminosity episode observed with 
INTEG RAL is similar to the on e s see n by Beppo- 



SAX dRea et al 


|2005[ 


& Galloway et alJ 


2000). 



such a low-luminosity episode showed a puls ed signal and 
a spec trum with an exponential roll-over. ICui fc Smithl 
l|2004h found with RXTE that, when the pulsation 
disappears or is very weak, the spectrum can be modelled 
by a simple power-law with index ~1.6. They also find 
episodes of suppressed flux that show a pulsation and 
a detectable roll-over in the spectrum. During the low- 
luminosity episode observed by INTEGRAL, the source is 
found to pulsate, but a simple power law with a slope of 
2.5 ± 0.2 models the spectrum of the source up to 150 keV 
without the need of an exponen tial cut-off. This seems in 
contrast with lCui fc Smithl l|2004h . On the other hand, we 
were not able to o bserve any signat ure of emission from 
a disk as found by iRea et all ( 2005h or an enha ncement 
of the absorption as found by Naik et ail l)2005(k this is 



due to the relatively high (17keV) lower threshold of our 



spectrum. In conclusion, with the analysed INTEGRAL 
observations we cannot settle the puzzle regarding the 
origin of these low-luminosity states. 

5. Conclusions 

In this paper we analysed the INTEGRAL data of the 
low-mass X-ray pulsar GX 1+4. The source was observed 
several times and found in different luminosity states with 
significant variations in its spectral and timing character- 
istics. 

The overall spin-down trend is confirmed by the high- 
est spin period ever observed for this source and by the 
negative local spin derivative we measured during most 
of the observation. However, in the high-luminosity state 
GX 1+4 showed a local spin-up. We estimate that the 
torque reversal happened at a sou rce luminosity simila r 
to analogous episodes described in IPereira et al.1 <jl 999ft . 
We notice that, during the spin up phase, the pulse pro- 
file is more peaked than in the intermediate luminosity 
state. 

When the source was brighter, the spectrum showed 
a sharper exponential cut-off after becoming harder be- 
tween 10 and 50keV. This is consistent with the thermal 
Comptonization being the dominant emission process. 

We also observed weak evidence of an absorption fea- 
ture in the phase-resolved spectrum o f GX 1+4 at 33- 
35keV in the phase interval [0.6-0.85]. foea et all <|2005h 
report a similar evidence in the phase-averaged spec- 
trum obtained with Beppo-SAX. If this feature is due to 
cyclotron resonant scattering, the pulsar magnetic field 
would be - 3 x 10 12 G. 
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Appendix A: ISGRI energy calibration issues 

The ISGRI detection layer consists of a 128x128 arra y 
of independent CdTe detector pixels ijLebrun et alJl2003]) . 
This kind of detectors is affected by the "Charge Loss 
Effect": since in CdTe the hole lifetime is much shorter 
than its transit time, the induced current can be signif- 
icantly reduced depending on the position of the inter- 
action on the detector. To reconstruct the energy of the 
incoming photons, it is necessary to perform an energy 
correction of the pulse amplitude value as a function of 
the pulse rise time; for this purpose OSA software makes 
use of multiplicative coefficients stored in a "Look-Up" 
Table (LUT2). 

The LUT2 distributed with the OSA software is not 
optimized, especially in an energy region near 80keV 
where a fraction of the photons is overcorrected. Ad hoc 
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absorption-like features at ~ 80 keV have been introduced 
in the OSA 4.2 and OSA 5.1 ISGRI Effective Area (re- 
spectively dashed and solid lines in Fig. IA. 1|) modelled 
on the basis of the observed Crab spectrum assumed to 
be a power law with index 2.2. This choice of the in- 
dex does not agree wi th the widely accepted value of 2.1 
l|Toor fc Sewardlfl974|) . ISGRI spectra are thus systemat- 
ically softer than other instruments. Moreover, the intro- 
duction of these ad hoc corrections in the effective area 
could produce the appearance of fake features in the spec- 
tra of sources that are different from the Crab. 

To improve the calibration, a new LUT2 based on on- 
ground and in-flight calibration data has been generated. 
Using this new LUT2 a better charge loss correction is 
obtained. 2 

We also generated new spectral-response matrices tak- 
ing into account the charge loss effect and its correc- 
tion by the new LUT2. The effective area, shown as the 
dash-dotted line in Fig lA.ll does not exhibit any 80keV 
absorption-like feature. This constitutes a confirmation 
of the better energy correction performed with the new 
LUT2: there is, in fact, no physical reason for this fea- 
ture. 

Also the spectral parameters we obtained with our re- 
sponse matrices are closer to the literature values (e.g. 
an index 2.1 for the Crab), while the x 2 values of the 
fits are systematically better than the ones obtained us- 
ing the OSA 4.2 and 5.1 calibration files. In Fig. lA~2l the fit 
to the Crab spectrum can be appreciated, but the reader 
can read a more d etailed discussion in the appendix of 
iMineo et all l)2006^ . 

If one assumes that the fractional spread of the Crab 
residuals of Fig. IA.2I is an empirical indicator of the cali- 
bration systematic errors, the calibration systematics can 
be neglected for our observation in OP4. In fact, the Crab 
reference spectrum we used and the source spectrum have 
about the same number of total counts in the 20-40 keV 
energy range but the S/N of the Crab is three times higher. 

Appendix B: Timing analysis of the ISGRI data. 

The standard OSA pipeline does not allow to perform 
timing analysis of pulsars. We then developed an alter- 
native software 3 to extract simultaneously spectra, light 
curves and phase resolved spectra from the event lists; 
a dedicated paper which describes the method in detail 
(Segreto et al., in prep.) is in preparation, here we give a 
brief overview. 

When illuminated by an X-ray source, a coded mask 
casts a shadow on the detector plane (shadowgram) that, 
correlated with the known mask pattern, allows to obtain 



OSA ARF 



ADOPTED ARF 



2 For the technical details we refer to Alberto Segreto's talk 
at the Internal INTEGRAL workshop available at 

http : //www. rssd. esa. int/Integral/workshops/Jan2005/ 
and to the new calibration files available at 
http : //www. if c . inaf . it/f errigno/ integral/. 

3 The software can be downloaded at the URL 

http : //www. if c . inaf . it/f errigno/INTEGRALsof tware .html. 
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Figure A. 1. Comparison between the effective area as 
distributed with OSA (left panel) and the new effective 
area obtained after performing the energy correction with 
the new LUT2 (right panel). In the left panel the dashed 
line is the effective area distributed with OSA 4.2, and 
the solid line the effective area distributed with OSA 5.1. 
The OSA effective areas present many bumps and wiggles, 
which are introduced ad hoc to compensate for the OSA 
LUT2 energy correction. Features in the 30-40 keV region 
are due to CdTe absorption edges. 
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Figure A. 2. Crab spectrum from a selection of science 
windows in the revolutions 102 and 103 in staring obser- 
vations, the exposure is about 47 ks. The fit is obtained 
with a power law (slope 2.1) without adding systematic 
errors. 



the position and in tensity of the sources w ithin the instru- 
ment field of view l(Coldwurm et al]l2003h . 

To obtain spectra and light curves of these sources, a 
shadowgram deconvolution process can be repeated sev- 
eral times, each time after selecting events in the energy 
and time intervals of interest. However, when narrow en- 
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ergy or time bins are considered, a large amount of CPU 
time is required. 

A much faster way to obtain light curves or spectra 
from coded mask telescopes when the source positions 
have already been determined (because known a priori 
or predetermined by a shadowgram deconvolution) , is via 
the use of the illumination pattern that the mask casts 
on the detection plane. Given the pixellated nature of the 
detector plane, some of the detector pixels will then be 
fully illuminated by a source while other pixels will be 
partially illuminated or not illuminated at all. For each 
pixel the amount of illumination from each source, nor- 
malized to the maximum illumination value is called the 
pixel Photon Illumination Fraction (PIF). 

Using this PIF information, spectra, light curves, or 
other scientific products can be extracted by means of 
a minimum amount of algebraic operations, without the 
need of repeating a shadowgram deconvolution process 
over and over again. This meth od, which is used e.g . 
by the standard pipeline of BAT JlBarthelmv et alJl2005j) 
on board Swift l|Gehrels et all l2004|) and bv the JEM- 
X software, does not account for the presence of multi- 
ple sources in the same field of view. Due to a non ideal 
imaging point-spread function this may produce cross con- 
tamination between the signals from different sources. To 
avoid this problem, we then developed a decontamina- 
tion algorithm, based agai n on simple algebraic operations 
l|Segreto k FerrignolEoO^) . 

To validate our extraction method we compared the 
background-subtracted spectra of several sources with the 
one obtained using the OSA standard pipeline and verified 
that the two methods give equivalent results not only in 
terms of spectral shape but also in their S/N. 

The PIF-based extraction method we developed is not 
only much faster, but it also allows to easily perform fast 
timing analysis and phase-resolved spectroscopy, tasks not 
supported by the standard OSA software. 

To estimate the reliability of the background rejec- 
tion we extracted spectra in several sky directions chosen 
around the main source to ensure the same observational 
conditions but far enou gh to avoid proble ms due to the 12' 
point-spread function l|Gros et al.ll2003|) . In Fig. IB. II we 
plot an example for five pointings within OP4. As shown, 
the spectra are all consistent with zero without points with 
S /N exceeding ~ 3 in each spectrum. We verified that the 
S/N distribution of the spectra in the 19-110keV range is 
described well by a Gaussian function with er ~ 1. 
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Figure B.l. Spectra extracted in five positions at two 
degrees from GX 1+4 for the OP4 data set. 
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